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Abstract 

We present cosmological solutions for (1 + 3 + n)-dimensional steady state universe in dilaton gravity with an 
arbitrary dilaton coupling constant w and exponential dilaton self-interaction potentials in the string frame. 
We focus particularly on the class in which the 3-space expands with a time varying deceleration parameter. 
We discuss the number of the internal dimensions and the value of the dilaton coupling constant to determine 
the cases that are consistent with the observed universe and the primordial nucleosynthesis. The 3-space 
starts with a decelerated expansion rate and evolves into accelerated expansion phase subject to the values 
of w and n, but ends with a Big Rip in all cases. We discuss the cosmological evolution in further detail for 
the cases w = 1 and w — \ that permit exact solutions. We also comment on how the universe would be 
conceived by an observer in four dimensions who is unaware of the internal dimensions and thinks that the 
conventional general relativity is valid at cosmological scales. 
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1 Introduction 

The accelerated expansion of the universe first came under scrutiny in 1917 with de Sitter's vacuum solution 
of the Einstein's field equations with a positive cosmological constant. However, it has not been considered 
seriously before the inflationary scenario [T] that is characterized by an epoch of exponential expansion in the 
very early universe at energy scales ~ 10 16 GeV was proposed in 1980. The inflation idea is one of the most 
prominent attempts to resolve the problems of standard Big Bang Cosmology such as the observed spatial 
homogeneity, isotropy, and flatness of the universe and others (See [5] for a recent review). The idea has 
many variations which are generally based on conventional general relativity theory together with a number 
of scalar fields introduced in an ad hoc way. Yet it is still a scenario as a concrete derivation of inflation 
from a fundamental theory such as the string theory [3] hasn't been found. On the other hand, the fact 
that our universe has entered a state of accelerated expansion at redshifts less than ~ 0.5 is well established 
with independent studies [H [7J |H1 HI EH HU H2]- The lack of a satisfactory physical explanation of 
the current acceleration is a mystery because it happens at energy scales ~ 10 _4 eV where we supposedly 
know physics very well. The observed dynamics of the universe could be accommodated surprisingly well in 
ACDM cosmology by the inclusion of a positive cosmological constant into the Einstein's field equations in 
the presence of pressure-less matter. But now the cosmological constant problem has to be faced which is one 
of the pressing unsolved problems in physics today [121 [T21 Q~5] : The observed energy density of the vacuum 
turns out to be non-zero but some 120 orders of magnitude less than the theoretical estimates. This is about 
the same order of magnitude as of the baryonic matter, a fact that makes the present epoch of the universe 
very special (coincidence problem). 

The efforts to solve cosmological problems in the context of conventional general relativity led many to 
consider various sources of dark energy. However, the known dark energy candidates such as quintessence, 
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k-essence, phantom, quintom and others are mostly ad hoc and/or considered phenomenologically rather than 
being derived from a fundamental theory [16 . Another approach to cope with the current acceleration of the 
universe is to consider the possibility that the Einstein's theory of gravity (viz. general relativity) is not an 
adequate theory for understanding gravitation at cosmological scales. There are many alternative theories of 
gravity suggested with extra fields (scalar-tensor, Einstein-aether, ...), including higher derivative and non- 
local terms in the action ( i? 2 , f(R), ...) or living in higher dimensions (Kaluza-Klein, Randall-Sundrum, 
...) (See |17L I18| for comprehensive reviews). In fact any successful unified theory of everything, such as a 
superstring theory, is expected to resolve all the unexplained issues once and for all. Moreover, cosmology may 
turn out to be the main avenue to probe string theories, unless supersymmetry is discovered at low energies 
and its properties provide some hints to its high energy origin or else the string scale turns out to be low 
enough to be probed at LHC or some future collider [3J. 

Motivated by the arguments above, we consider here a higher dimensional dilaton gravity theory that 
may arise as the low energy effective field theory approximation to a superstring theory. It is well known that 
bosonic string theories require a 26-dimensional space-time for consistency while anomaly-free superstring 
theories live in 10 dimensions. It is generally assumed in the higher dimensional context that all but four 
of the space-time dimensions are compactified on an unobservable internal manifold, leaving an observable 
effective action to be considered in a (1 + 3)-dimensional space-time [191120) . There is yet another approach to 
the problem of extra dimensions. They may be relaxed to become dynamical so that they start at larger scales 
in the early universe and evolve in time to unobservable scales, thus leading to an effective four dimensional 
picture of the universe (See [21j for a comprehensive review on Kaluza-Klein cosmologies). This approach 
has been considered for the first time in the early 1980s as a dynamical reduction of internal dimensions with 
the external dimensions expanding while the internal dimensions contracting |22 l I23L [2"1] . Among the higher 
dimensional cosmological models, there is an interesting subclass that describe higher dimensional steady state 
cosmologies with the total volume of the higher dimensional space remaining constant in time [231 
|2"T] . In these papers, the external space expands with a constant deceleration parameter which implies either 
a power- law expansion or exponential expansion. In a recent study [28J , we have assumed from the beginning 
that both the higher dimensional volume of the universe and the energy density are constants and found 
some new solutions where the external space exhibits dynamical behavior with time dependent deceleration. 
We further concluded that the model may be improved by considering modified gravity theories. Indeed, 
in this paper considering a higher dimensional dilaton gravity we not only present power-law dynamics but 
also dynamics in which the external space accommodates the observed universe starting from the primordial 
nucleosynthesis epoch to the present universe. In all the models that can accommodate the observed universe, 
on the other hand, the external space ends with a Big Rip and hence lives for a finite time. This is also 
a result consistent with the observational studies that argue the Big Rip behavior for the future of the 
universe is slightly favored when compared, for instance, with the de Sitter future of the ACDM cosmology 
pM I3TJ1 I3T1 13"2l 155] , Moreover, because the universe lives for a finite time the coincidence problem is not a 
severe problem anymore as discussed first by Caldwell et al. |34| . 

The organization of the paper is as follows. In Section 2, we introduce the dilaton gravity action in 
(1 + 3 + ri)-dimensions written in the string frame. The variations of the action and the cosmological field 
equations are also given. We discuss in Section 3 steady state cosmological models with a constant (3 + n)- 
dimensional volume. We give two classes of exact solutions with the external space expanding while the 
internal space is contracting in Section 4 and Section 5 and study their observational consequences separately 
in detail. A brief discussion on how the universe appears to an observer living in 3-dimensions is given in 
Section 6. We conclude by some final comments. 



2 



2 Dilaton Gravity Equations 



Inspired by the low-energy effective string theory, we consider the following dilaton gravity action in (1 + 3+n)- 
dimensions written in the string frame: 

1=1 d 1+3+n x^\e- 2tfi {TZ + 4w dnipWip + U(tp)} . (1) 
J M 

1Z is the scalar curvature of the space-time metric g and tp is the dilaton field. The dilaton self-interaction 
potential U(ip) is a real, smooth function of the dilaton to be specified. We introduce an arbitrary positive 
dilaton coupling constant w > to cover a variety of cases in the string frame. For instance choosing w = 1, 
we obtain gravi-dilaton string effective action with a dilaton self-interaction potential in (l + 3 + n)-dimensions 
and the further choices 1 + 3 + n = 10 and 1 + 3 + n = 26 correspond to anomaly-free superstring and bosonic 
string theories, respectively. 

We consider a spatially homogeneous but not necessarily isotropic (1 + 3 + n)-dimensional synchronous 
space-time metric that involves a maximally symmetric three dimensional flat external space metric and a 
compact n-dimensional flat internal space metric: 

ds 2 = -dt 2 + a 2 (t) (dx 2 + dy 2 + dz 2 ) + s 2 (t) (d<9 2 + ... + d6> 2 ) . (2) 

Here t is the cosmic time variable, (x, y, z) are the Cartesian coordinates of the 3-dimensional flat space that 
represents the space we observe today, and (9i,...,0 n ) are the coordinates of the n-dimensional, compact 
(toroidal) internal space that represents the space that cannot be observed directly and locally today. The 
scale factors a(t) of the 3-dimensional external space and s(t) of the n-dimensional internal space are functions 
of t only. 

Here we will be considering an exponential potential for the dilaton field 

U(ip) = C/ e V , (3) 

where Uq and A are two real parameters and the dilaton field is space independent real function of time t: 

<p=<p(t). (4) 

The field equations to be solved are determined from the action ((TJ) by a variational principle. The usual 
way to proceed at this stage is to substitute the cosmological field expressions into these field equations and 
reduce them to a system of coupled ordinary differential equations. Here we will proceed in the following way 
to reach this final set of equations to be solved. Let us first substitute the cosmological fields into the action 
integral so that the reduced action reads 

I[a, s,<p] = J dt J d 3 x J d n 9C (5) 

where the Lagrangian function 

£ = e- 2v a 3 s n {n- iwtp' 2 + U e Xv } (6) 

with the curvature scalar given by 



a'V a" , , ( s'\ 2 _ s" . a' tt 



U = 6 — + 6— + n(n - 1) — + 2n— + 6n 

a J a V s / s a s 

Here ' denotes 4z. Differentiating by parts and neglecting total derivatives we arrive at an equivalent La- 
grangian function 

C = e- 2 ^a 3 s n \-6n--- 6 ( -) - n(n - 1) ( -) + + 4n-cp' - 4w<p' 2 + U e^ \ (7) 

a s \ a J \ s J a s 
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from which we get the following Euler-Lagrange equations of motion 

n' 2 n" n' s' /I 3 \ s' 2 s" 

3^ + 3— + (3n - 3)-- + V-^ + lU-i^-l)-- 2^" (8) 

+(4 - 2uV 2 - 2^ (^3- + (n - + = 0, 

\ a s J 2 

a' 2 a' s' n(n-l)s' 2 a" s" „ , 2 , / a' s'\ 2- A Am 

3^- + 3n + — ^- — ~ + 3— + n 2wip" + 2ww' - 2ww' 3— +n— + — ^^oe v = 0, (9) 

as 2 s z a s ' \ a s J 4 

a' 2 a" oV s" n(n~l)s' 2 ,. , s , 2 , / a' s'\ 1 Am 

2y>" + (4 - 2w>p' -2</j' 2— +n- ) + -U e Xip = 0, (10) 



., , 2— + 2n +n— , , _ r _ . , r _ r 

a^aass 2 s z \ a s / 2 

for the dynamical variables a, and s, respectively. These coupled equations must be supplemented by the 
Hamiltonian constraint equation 

3^ + in*-*- + fcll ^ + 2V 2 _ V ( 3 ^ + n S l) + = . (11) 

a z as 2 s z \ a s J 2 

We would like to note that any one of the equations out of (|5])- (fTTjl can be derived from the remaining two 
other equations plus the Hamiltonian constraint equation (See, for instance, |35| for details). 



3 Cosmological Models with (3 + n)-Dimensional Constant Volume 

We are interested in the solutions for a higher dimensional steady state universe. In a recent study [28], we 
define a higher dimensional steady-state universe with two properties: (i) the higher dimensional universe has 
a constant volume as a whole but the internal and external spaces are dynamical, (ii) the energy density is 
constant in the higher dimensional universe. However, we do not introduce a matter source here and hence 
it is enough if we consider only the first one of the properties. Accordingly, we assume that the volume scale 
factor of the higher dimensional universe, i.e. the (3 + n)-dimensional volume, is constant: 

V = a 3 s n = V . (12) 

Vq is a positive real constant. We note that this condition assures the dynamical contraction (hence reduction) 
of the internal space for an expanding external space and that the constant (3 + n)-dimensional volume may 
be interpreted either as the total volume of the (3 + n)-dimensional space of the (1 + 3 + n)-dimensional 
universe or as the volume of a (3 + n)-dimensional space-like domain in a (1 + 3 + n) dimensional universe. 
Substituting (3 -I- n)-dimensional constant volume constraint (fT2"j) in the equations (j5|- (fTTj) , they reduce to 
the following set of equations: 

|f- + l)4 + --- V + (4 - 2^' 2 - V- + = 0, (13) 

2 \n J a z n a n a 2 

2 wtp" + 2 V 2 + 2 —^U e^ = 0, (14) 
2 \n J a 1 4 

IT" + I) 4 - - - V + (4 - 2uV 2 + 2^' + \u,e^ = 0, (15) 
\ 2n 2 J a z a a 2 

-U- + l)^ + lW 2 + \u^ = Q. (16) 
2 \n J a 1 2 

Now adding equations (|T4l and (TTB1 side by side, the dilaton field decouples from the scale factors and satisfies 

16?V 2 + (4 - X)U e Xip - Swtp" = 0. (17) 
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Furthermore subtracting (fT5|) from fT3")) and then using (fTT)) we obtain 

, la" la' 
2 a 2 a 

which finally gives 

a = a Q e Cl f e2 * dt (19) 

where ao and c\ are real integration constants. Now subtracting (|14j) from (|15p and substituting (If 91) we 
obtain 

cm(A - 4 + Aw) (2ip' 2 - y"\ = 0. (20) 

We are looking for solutions where the universe is possibly dynamical. Then c\ ^ from (|19p and we must 
set n 7^ 0, otherwise there will no internal dimensions and a 3 = Vb, means that the universe is constrained 
to be static. Hence a dynamical cosmological model with constant (3 + n)-dimensional volume scale factor 
constraint is not consistent for arbitrary values of the constants w and A, but under the condition 

\ = A-Aw. (21) 

This condition reduces (|17fl to the following second-order non-linear ordinary differential equation 

if/' - 2<p' 2 - = o. (22) 

Now defining a new scalar field 

/3(i)= e - 2 ^, (23) 

which is a positive function since <p is a real function, and using this in (|22p we obtain the following second- 
order autonomous differential equation 

/3" = -UqP 2 ™- 1 . (24) 

Integrating this equation once we obtain 

2£ = P' 2 + — B 2w (25) 
w 

where E is a real constant of integration. This differential equation is familiar from the classical central force 
problem in one dimension in the presence of a power-law potential, though with a slight difference. One may 
observe from (|2"3"]) that j3 is non- negative j3 > 0, since dilaton field cp is defined to be real function. This should 
be kept in mind in what follows as we proceed with the solutions considering /3 > 0. We would like to expose 
the meaning of the integration constant £. Equation (|25p can be rewritten as 

£ = \f - J -Uap 2 ™- 1 ^. (26) 

It is of the form 

S = T - J F{P)&(3 = T + U(fi), (27) 
where T is the kinetic energy and U is the potential energy given by 

T = -P' 2 and U(B) = - f HPW = —P 2w , where T{B) = ~~U a p 2w -\ (28) 
2 J 2w 

Hence, £ corresponds to the total energy of the j3 field. From (|23|1 we have 

d/3 



28- LIslr^ 



= ±t-t , (29) 
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where to is an integration constant and can be chosen arbitrarily since it is nothing but a shift in the cosmic 
time t. For convenience we choose to = and positive sign for the time parameter t. 

Finally, a relation between the scale factor of the external space a and the f3 field (hence dilaton field ip) , 
including the dilaton coupling constant w, the total energy £ and the number of the internal dimensions n 
can now be given: 

^ = «£ttw-- (30) 

a z 3(3 + n) p z 

We note that £ > 0, since w > 0. Integrating ((3T)|) we obtain the scale factor of the external space, and hence 
of the internal space via (|T2"j) . as follows 



a = aoe ±i^^f^t and s = Soe ^^^f^\ (31) 

oo and so are integration constants that satisfy constant volume condition Vq — o>q s o- Hence, the scale factors 
may be obtained explicitly whenever /3 can be obtained explicitly from (j29[) . However, in general it is not 
possible to evaluate the integral in (|29[) in terms of known functions for arbitrary values of w, Uo and £ , but 
only for particular choices. It is possible to evaluate the integrals and isolate /3 for arbitrary values of w in two 
cases: (i) when the total energy of the (3 field is zero £ — and the potential of the dilaton field is non-zero 
Uo 7^ 0, (ii) when the dilaton field has no potential Uo = and the total energy of /3 field is positive £ > 0. 
On the other hand, when £ > and Uq ^ 0, it is possible to evaluate integral and isolate j3 for w = 1, i; for 
other values of w > it is possible to obtain results for /3 in terms of circular or elliptic functions but in many 
cases only numerical results maybe gained. Fortunately, we are able to get some information for arbitrary 
values of w > that may be important from the cosmological point of view. To do so, we first give some 
important cosmological parameters in the terms of j3 and its derivatives with respect to time. Using pip we 
obtain the Hubble parameters of the external and internal spaces, respectively, as follows: 



a' 1 / 3n V2w£ , „ s' 1 / 3n V2w£ 
H a =-=±-J V and H.= - = T-J — V . 32 

a 3 V 3 + n p s n V 3 + n p 

We note that Hubble parameters are inversely proportional with j3 and either the external space expands 
while the internal space contracts or conversely the external space contracts while the internal space expands 
such that (3 + n)-dimensional volume scale factor is constant, i.e., Vq/Vq = 3H a + nH s = 0. Because we 
are interested in the solutions where the external space expands in particular in what follows we shall dwell 
on solutions with expanding external space, i.e., the cases with H a > 0. The dimensionless deceleration 
parameter is another important cosmological parameter and the deceleration parameters of the external and 
of the internal dimensions are 



a"a ^ /3 + n p s"s 3 + n ? 

q a = =- = — 1±S\ , and q s = — = —i=pn\ , (33) 

a' 2 V 3n y^£ s' 2 V 3n V2^£ 

respectively. The deceleration parameter is particularly important when we are discussing the dynamics of 
the external space that represents the space we observe today. The negative values of deceleration parameter 
correspond to the acceleration and positive values to the deceleration of the scale factor. The space is said 
to exhibit super-exponential expansion for q < — 1; exponential expansion (also known as the de Sitter 
expansion) for q = — 1; accelerated power- law expansion for — 1 < q < 0; expansion at a constant rate for 
q = 1 and finally decelerated expansion for q > 1. In the ACDM cosmology, which is the simplest model 
that accommodates the observed dynamics of the universe but with some serious problems, the deceleration 
parameter evolves from i to —1 as the universe expands |15| . In most of the scalar field dark energy models 
in the context of (1 + 3)-dimensional general relativity, on the other hand, it evolves between 2 and —1 |16) . 
Because we are concerned in this paper with an evolving deceleration parameter, it is convenient to discuss 
also the dimensionless jerk/state finder parameter (that involves the third derivative of the scale factor of the 
observed universe with respect to time) of the external space into account: 



,3 + n ( „ „, „ /3 + n 
(/' a \ v in 



j a ^-^= ±6^ 2/3' - 9\I^Jlf3> - + l)w£. (34) 
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The jerk parameter in ACDM is simply equal to unity j'acdm = 1 |36J. Hence the search for variation of j from 
unity over some redshift interval or time interval is also very convenient for seeing the departures from ACDM 
model. However, the observational constraints on the value of jerk parameter are rather weak in the present 
literature. The various observational studies give the observed value of the jerk parameter approximately in 
the range -5 < j , to day < 10 [3H EH EH EH1 SHI US S21 H3] ■ 



3.1 The Static Universe 

The first case we would like to discuss is the one where £ = and Uq ^ 0, i.e. the total energy density of the 
(3 field is null £ = but the potential is non-zero Uq ^ 0. However, we note from (|25p that the case £ = is 
possible only if Uq < 0. Now choosing £ = in (|2T)|) we obtain: 



(1 - wfU Q ^ 



w 



and ip 



4(1 - w) 



In 



(1 - w) 2 U f2 



w 



and 



P 



' Uot and <p=}Q?°t for w = 1. 



for w > 1 



(35) 



(36) 



On the other hand, although the dilaton field is dynamical, the universe is non-dynamical in this case. 
Choosing £ — in (l3"Tjl we have 

a = ao and s = sq. (37) 

This result tells us that if £ = 0, that is, if the total energy of the f3 field is zero, then we have a static universe 
independent of the value of dilaton coupling constant w. 



3.2 Dynamical universes 

Before giving the solutions with dynamical external and internal spaces we would like to give a short discussion 
on whether we can have consistent effective 4-dimensional cosmologies or not under the assumption of constant 
higher dimensional volume independent of a specific solution. 

We first comment on the role of the number of internal dimensions n on the expansion rate of the 
external dimensions and the contraction rate of the internal dimensions and their sizes. From the constant 
(3 + n)-dimensional volume condition, the magnitude of the expansion rate of the external space and of the 
contraction rate of the internal space are identical independent of the number of the internal dimensions. 
However, one may observe from (|32p that the number of the internal dimensions n appears as two different 
factors in the Hubble parameters of the external and internal dimensions and affects the expansion rate of the 
external dimensions slightly but the contraction rate of the internal dimensions drastically so that for very 
large n values the internal dimensions are almost frozen. Apparently, this has consequences on the size of 
the external and internal dimensions. From the constant higher dimensional volume condition, the size of the 
internal dimensions can be given as follows: s — Vo n a _ «. This relation gives us opportunity to investigate 
whether the internal dimensions are involved with the process of primordial nucleosynthesis, in the early 
universe, that shouldn't be affected by the presence of internal dimensions for any viable cosmological model. 
To do so, we assume that the size of the internal dimensions remain roughly less than the size of a proton 

^proton ^ 

1 GeV 1 ~ 10 16 m when the primordial nucleosynthesis occurs, i.e., sns ^ proton- Additionally we 

can write 

S =(^)" Stoday (38) 

from the constant higher dimensional volume condition (fT2"j) . As long as the wavelength of the cosmic 
background radiation (CBR) is so high that it cannot propagate into the internal dimensions s (m) < 
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438 T- 1 (Kelvin" 1 ) we can use the standard relation^ between the CBR temperature and the scale factor 
of the external space 

^today J- 



0* ^today 

Accordingly a successful model should provide the condition 



(39) 



S NS = ( „ NS ) Stoday < 'proton ~ 10 16 m. (40) 
\ J today / 

The primordial nucleosynthesis of the light elements is determined by the events occurring in the epoch 
when the temperature of the universe varied from ~ lMeV (~ 10 9 K) to ~ 0.1 MeV (~ 10 s K) and the 
current temperature of the CBR is T tod ay = 2.352 x 10~ 4 eV (T tod ay = 2.728 K) [ill [12], Accordingly, we 
find that atoday/«NS ~ 10 9 from and Stoday/^NS ~ 10~~ from (|4"0|) . We first check if the constant 
higher dimensional volume constraint remains valid under an extremely strong condition such as the case 
when the size of the internal dimensions today are almost at LHC length scales s to day ~ 'lhc ~ 10 _20 m 
('lhc ~ 10 _13 eV _1 ) and can be observed very soon. In this case the condition sns ^ 'proton is satisfied 
provided n > 7. On the other hand, assuming the size of the internal dimensions today are at Planck length 
scale Stoday ~ 'pianck ~ 10 _35 m ('pianck ~ 10 _28 eV _1 ) then sns ^$ 'proton condition is satisfied provided n > 2. 
We can also estimate the mean length scale of the volume of the higher dimensional universe I = (a 3 s n ) 3+™ 
using sns J$ 'proton for Stoday ~ 'lhc and Stoday ~ 'pianck- The size of the observed universe today can be 
taken as a to day ~ 10 26 m. Using Stoday ~ 'lhc and n > 7 we obtain I < 10~ 6 m and using Stoday ~ 'pianck 
and n > 2 we obtain I < 10° m for the mean length scale of the higher dimensional constant volume. The 
analysis above shows that we can have cosmological models with higher dimensional constant volume without 
the intervention of the internal dimensions to the primordial nucleosynthesis processes. Moreover if we can 
obtain the scale factor of the external space a explicitly, we will be able to examine the physical processes 
in the early universe in detail using (|39p . Hence one can further investigate if the dynamics of the 3-space 
allows a primordial nucleosynthesis process in the early universe (see |45l |4"tj] for details). To ensure that 
some primordial nucleosynthesis in the early universe takes place the temperature of the CBR should reach 
T ~ 80keV before the free neutrons decay t < 887 s. The standard Big Bang nucleosynthesis occurs when 
the temperature varies from T ~ 1 MeV to T ~ 0.1 MeV, meanwhile the universe evolves very similar to the 
power-law cosmology as a ~ £~o.55 which gives a deceleration parameter q ~ 0.82 [46]. Hence we are able to 
check whether a model allows primordial nucleosynthesis and further if it produces the light elements with 
correct abundances using (IBTJt once the scale factor of the external space is obtained explicitly. 



3.2.1 The case 8 ^ and U = 

This is the pure dilaton solution since the dilaton potential is null Uq = 0. We note from (|25|) that a solution 
is possible only if £ > 0, and now choosing Uq = in ([29)) we obtain 

(3 = V2£t and y = — In (V2S tj , t>0. (41) 

The scale factors of the external and internal dimensions are as follows: 

a = a t ± iV^^ and s = s t T - v^v^ (42) 

This is the well-known power-law dynamics and yields the following constant deceleration parameters for the 
external and internal dimensions: 

l3 + n 1 , /3 + n 1 , . 

q a = ±3W — = - 1 and q s = ^pn\ — = - 1. (43) 

V 3n y/w V on ^Jw 

In fact the relation is as follows Qtoday = Ptt-^- — , where /3 accounts for non-adiabatic expansion due to entropy 

a -'today 

production. In standard cosmology ft = 1 for T < m c i octr on and j3 — (t) s f° r T > "icicctron- However, we simply use 
P = I in this study, since it is a good approximation within the context of our discussions. See Ref. |46| and references 
therein for details. 
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In this solution the expansion/contraction of the external space and contraction/expansion of the internal 
space start at t = and last forever. We note that the values of the deceleration parameters are determined 
together by the number of the internal dimensions n and dilaton coupling parameter w and accelerated 
expansion (q a < 0) is also possible for the external space. One may observe that the external space exhibits 
accelerated expansion, if w > 12 for n = 1, and if w > 3 in the n — > oo limit. On the other hand, 
accelerated expansion of the external space is not possible if w < 3, independent of the number of the internal 
dimensions. To be precise considering the present value of the deceleration parameter of the observed universe 
Qa, today = —0.81 we find that 



As the presence of a matter source does not change the dynamics considerably and the internal dimensions 
are small enough, we find that the models for 



can produce correct abundances of the light elements. Accordingly, as we expected, the power-law dynamics 
cannot represent the complete picture of the universe we observe but may be valid during a particular period 
of the universe depending on w and n. Moreover, the cost of representing the currently accelerated universe 
is very high w values. Then, the question we face is whether solutions that give a more complete picture of 
the observed universe for a specific value of w value, and hence solutions where the external space exhibits 
a variable deceleration parameter consistent with the observed universe can be obtained. We give promising 
solutions in this respect in the following sections. 

3.2.2 The case £ ^ and U / 

This is the case where we obtain a time- varying deceleration parameter and hence is the one we are particularly 
interested in. First of all, we note from ([3H)l that solutions exist only if £ > as we have already discussed 
above. Then one may observe from ([2U)) and (fBTfl) that we have various dynamics depending on the value of the 
dilaton coupling constant w since it determines the behavior of the j3 field. We are particularly interested in 
the solution where the observed universe (the 3-dimensional external space) is expanding and its deceleration 
parameter is a decreasing function of time, i.e. q' a < 0. One may observe from (f3"3")l that q' a < is satisfied 
provided that f3" < 0. This condition also implies from that Uq > 0. Hence in what follows we discuss 
solutions for £ > and Uo > only. 

We observe from (I5T|) that the scale factor of the expanding external space goes either to zero or to 
infinity depending on the sign of the exponent as /3 — > 0. It follows from (|25p that /3 = at /3' = ±\/2£. Then 
rewriting the integral as J ^ = J d ^" ^ , we see that it goes to minus infinity (hence a — > 0) as /?' — > \p2£ and 
goes to plus infinity (hence a — > oo) as /?' — >• —\f2£ independent of the dilaton coupling constant w. On the 
other hand, the evolution trajectory of the function j3' and hence of the scale factors are determined by the 
dilaton coupling parameter w. However, note that the value of f3' at f3 = is determined only by the value 
of £ . Moreover, f3' is the only function that appears in the deceleration parameters ([3"3")) . Therefore we are 
able to explore the initial and final values of the deceleration parameter of the external dimensions, although 
we cannot obtain the function explicitly for arbitrary w values. Thus the deceleration parameter of the 
external dimensions when a = (/3' = V2£) and a — > oo (/?' = — \[2£ ) from are as follows: 




(44) 




(45) 



ensure some primordial nucleosynthesis, and that the models for 




(46) 




and 




(47) 
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We note that the potential of the /? field is null U = when j3' — \/28. Hence q a =o we obtained here is the 
same with the one we obtained for the case £ > and U$ = 0, i.e. the case with zero potential, in the previous 
section. Hence the initial (q a =o) and the final (q a =oc) values of the deceleration parameter of the external 
space are determined by the number of the internal dimensions n and the value of the dilaton parameter w 
only. One may observe that q a =o > — 1 and q a =oo < — 1, and hence the evolution of the external space passes 
through q a = —1 (when f3 reaches its maximum, i.e., j3' = 0) for any finite value of w > 0. 

Because of the fact that we observe a universe which is evolving from a decelerated expansion phase to an 
accelerated expansion phase we expect at the beginning the deceleration parameter of the external space to 
be larger than at a — 0. Therefore, we are able to put limits on the value of w with respect to the number 
of internal dimensions n: 

q a =o >0 < w < 3 + -. (48) 

n 

Thus given the number of internal dimensions, we have a certain range of allowed w values that satisfy 
q a =o > ( see Fig- [IJ- The higher the number of internal dimensions the lower the maximum value of w. 
Hence choosing the lowest value for the number of the internal dimensions n — 1 we see that a 3-space that 
evolves from decelerated expansion to accelerated expansion is possible only if w < 12. For example, if n = 6 
then only the models with < w < | give a 3-space that starts expanding with a decelerated expansion rate. 
This inequality, in return, also puts constraints on the number of internal dimensions n. For example, given 
w = 6, the number of the internal dimensions should be less than 3. On the other hand, for w = 1 the 3-space 
starts expanding with a decelerated expansion rate for any number of the internal dimensions. 

We can tighten the bounds on the parameters by introducing upper and lower limits on q a =o so as to have 
an idea whether a solution can accommodate the observed universe. In the conventional four dimensional 
cosmology in general relativity, stiff fluid was suggested as the most promising equation of state (EoS) of 
matter at ultra-high densities for representing the very early universe [47l 148) . Stiff fluid can be represented 
with an EoS p — p, where p and p are the energy density and pressure respectively. Such an EoS gives a 
deceleration parameter equal to 2 in general relativity. Motivated by this observation we choose q a =o — 2 as 
the upper limit for a possibly viable cosmological model. For the lower limit we can simply choose q a =a = ^ , 
considering ACDM which is the simplest model that is consistent with the observations and in which universe 
evolves from q a =o — \ to <Za=oc = — 1- Hence the possibly viable values of n and w should satisfy the bounds 

1 „ 11 4 4 

^<<Za=o<2 =>■ - + - <w<- + -, (49) 

2 6 n 6 n 

which is also depicted in Fig. [5] For instance, if n ~ 6 possibly viable models should satisfy i < w < 2. On 
the other hand, if w = 2 only the models for n < 6 are viable or if w = 1 only the models for n > 2 are 
viable, w = | is the only value of the dilaton coupling constant that is allowed for any number of the internal 
dimensions, as giving q a= o = 2 for n = 1 and q a =o = \ for n — > oo. 

Checking (|49l) (or see Fig. \2§ we conclude that one may have models that describe the observed dynamics 
of the universe successfully if 0.33 < w < 5.33. Moreover, one would like the models to accommodate also 
the primordial nucleosynthesis processes. We note that this solution with dynamical deceleration parameter 
approximates the power-law solutions when a ~ hence in the very early times, say at times of the primordial 
nucleosynthesis. Therefore considering the bound we determined in the previous section for the power-law 
solution we observe that the models in the range 0.33 < w < 4.04 may allow primordial nucleosynthesis 
processes and that the models in the range 0.91 < w < 3.64 may allow primordial nucleosynthesis processes 
that further yields acceptable light element abundances. In short we have various models according to the 
value of the dilaton coupling constant w and the number of the internal dimensions n and can even have 
cosmological models in which effectively four dimensional universe that can provide suitable conditions for a 
successful primordial nucleosynthesis in the early universe and evolves from a decelerated expansion phase to 
an accelerated expansion phase consistently with the data from cosmic microwave background and supernova. 

We would also like to comment on the future of the universe in these models. It is well-known that the 
deceleration parameter values lower than —1, i.e. expansion rates faster than the exponential expansion, cause 
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Figure 1: The initial (q a -o) and final (q a =oo) values of 
the deceleration parameter of the external space versus 
dilaton coupling parameter w for different number of in- 
ternal dimensions n. We plotted for n — 1, 2, 6, 22 and for 
n — ¥ oo limit, in the order from the solid blue curve to the 
dotted red curve. 



Figure 2: The initial (q a =o) value of the deceleration 
parameter of the external space versus dilaton coupling 
parameter w for different number of internal dimensions n 
for 0.5 < q a =o < 2. We plotted for n = 1,2,6,22 and for 
7i — 5- oo limit, in the order from the solid blue curve to the 
dotted red curve. 



a scale factor to expand to infinitely large values in finite time, a behavior known as the Big Rip [34 . Hence 
in the models given in this section the external space will expand to infinitely large values in finite time since 
<7a=oo < — 1 for any combination of n and w. The same result can also be derived from (|32j) . H a starts with 
infinitely large values when a = (/? = and (3' = \[2E') 1 decreases and reaches its minimum when (3 = /3 max 
and starts to increase and becomes infinitely large as a — > oo (/3 = and (3' = —^/2£). 

Up to this point we have discussed the cases for £ > and Uo > qualitatively because j3 function is, 
in general, cannot be expressed in terms of known functions. However, it is expressible in terms of known 
functions and J \&t is integrable so as to obtain the scale factors explicitly in two cases: 

• Case w = 1 (A = 0): We let w = 1 in flU} so that 

0' = -U o 0, (50) 
which is the familiar simple harmonic oscillator equation and the solution of (|29[) gives 

13 = J— sm(^lTot), 0<t<^=. (51) 

• Case w = \ (A = 2): We let w = | in (|24]) so that 

P" = -U , (52) 

which is similar to the equation of motion of a particle under a constant attractive force in classical 
physics, and the solution of (|2lH gives 



B- £ l U V2£ V2£ 



In what follows, we discuss each of these two exact solutions in further detail. 
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4 The model for w = 1 (A = 0) 



The special choice w = 1 reduces the action ([T]) to the gravi-dilaton string effective action in (1 + 3 + n)- 
dimensions and explicit solution will be provided £ > and Uq > as discussed above. Now, using (|51[) in 
((23)) we obtain the dilaton field as follows: 

9 = -^hiU^sm(^t)\ 0<i<-^=. (54) 

Because A = for w = 1, the dilaton field yields a constant self-interaction potential: 

U = U . (55) 
The scalar curvature of the higher dimensional spacetime: 

K = 5 7° r . (56) 
sin 2 (v^ t) 

We obtain the scale factors of the external and internal as 

a = ai6 3\/^arctanh(-cos(v^t)) and g = Sl g-iy^arctanh(- cm(v^o *)) ) (57) 

respectively, where ai and si are integration constants that obey ai 3 si n = Vq. The Hubble parameters of the 
external and internal spaces are 



g 1/3?L f*_, and g a = -l./J!L (5g) 
3 V 3 + n sin( v / Do *) n V 3 + n sin(v^ *) 

The deceleration parameters of the external and internal spaces are 



/ 3 _|_ n / 3 + n 

go = -1 + 3W 3 ^ cos(V Uq t) and q a = —l — ny — cos(y/U t). (59) 

We also give the jerk parameter for the external space only, since it is important while discussing the model 
from the observational point of view: 



3a = ( 3+- )cos 2 (VC^ot)+9- l /^i^cos(V^o<) + -+4. (60) 
\ n J V on n 

The evolution of the universe starts with a zero size external space and infinitely large internal space at 
t = 0. The external space that we observe today expands until it becomes infinitely large at a finite time 
f end = -S=, while the internal space contracts to unobservable scales and then reaches zero size at t e nd- 

The external space enters into the accelerated expansion phase at t acc = arccos (h^J 3+^7) ■ passes over 
into the super-exponential expansion phase at t sc — 2 J V and then reaches infinitely large expansion rate at 
t = i C nd- We note that all the cosmological parameters are exactly determined by the values of Uq and n 
only; in particular Uo determines the lifetime of the universe, n determines the initial and final values of the 
deceleration parameters and they together determine the time of the onset of the acceleration of the external 
space. Hence, in principle, we are able to obtain the number of internal dimensions n and Uo upon obtaining the 
values of the deceleration parameter at two different times or cosmic redshifts from observations. However, 
in this paper, we won't undertake an observational study but discuss the dynamics for different number 
of internal dimensions simply by considering the present age of the universe and the present value of the 
deceleration parameter that are consistent with various observational studies, in particular concerning model 
independent studies. Hence, with the choice q afi = -0.81 [371 ES 0M \M El 111 \M and ^today = 13.7 Gyr pj] 
considering the latest observations, we are able to depict the dynamics of the universe for different number of 
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Table 1: Cosmological parameters for the case w = 1 (A = 0) for different number of internal dimensions n assuming 
q a , o = -0.81 at ttoday = 13.7 Gyr. 





n = 1 


n = 2 


n = 6 


n = 22 


n — > oo 


•2 acc 


0.07 


0.12 


0.22 


0.31 


0.37 


13.7 Gyr - i acc 


2.15 Gyr 


2.78 Gyr 


3.71 Gyr 


4.38 Gyr 


4.74 Gyr 


^end 


28.39 Gyr 


28.67 Gyr 


29.06 Gyr 


29.32 Gyr 


29.46 Gyr 


Qa,int-> Ja,int 


2.46, 14.61 


1.74, 7.78 


1.12, 3.64 


0.85, 2.28 


0.73, 1.80 


Qa,0: ja,0 


-0.81, 12.47 


-0.81, 7.97 


-0.81, 4.97 


-0.81, 3.88 


-0.81, 3.47 


Qa,cndi Ja, end 


-4.46, 35.39 


-3.74, 24.22 


-3.12, 16.36 


-2.85, 13.36 


-2.73, 12.20 


so 


< 10- 43 m 


< 10" 30 m 


< 10~ 21 m 


< nr 18 m 


< 10" 16 m 




10" 12 s < 887 s 


10~ 6 s < 887 s 


1Q- X s < 887 s 


94s < 887s 


884 s < 887 s 



internal dimensions n. Accordingly, we plot the scale factors of the external and internal dimensions in Fig. [3] 
and the Hubble parameters in Fig. 2]for n = 1, 2, 6, 22 and also for n — >• oo limit. We plot the dimensionless 
deceleration parameter of the external dimensions in Fig. [5] and the jerk parameter in Fig. [5] The evolution 
of the deceleration parameter with the cosmic redshift is important to have an idea about the evolution of 
the observed universe. Hence, we also give the deceleration parameter of the external space versus the cosmic 
redshift z = — 1 + (where a z= o is the present value of the scale factor of the external space) 

q a (z) = -1 + 3^±Ztanh (s^m^ ln(l + z) + arctanh U^^-(q + 1)J ^ (61) 

and plot the deceleration parameter of the external space in terms of cosmic redshift z in Fig. [7] We 
summarize, in Table [JJ the plots by giving the cosmic redshift (z acc ) and the time that passed since the time 
of the onset of acceleration (13.7 Gyr— i acc ), the time of the Big Rip (tend) an d the initial (q a ,int , Ja.mt), present 
(q a .o, jafl) an d final {q a ,endi ja,ead) values of the deceleration and jerk parameters. In Table[JJ assuming that 
the presence of a matter source does not alter the dynamics considerably, we also calculated the time tx^so keV 
to see whether T ~ 80 keV < 887 s is provided to ensure some primordial nucleosynthesis. The maximum size 
of the internal dimensions today sq is also given considering the condition sns ^$ ^proton in Table QJ 

We observe that while the number of the internal dimensions alters the expansion rate of the internal 
dimensions slightly, it alters the contraction rate of the internal dimensions drastically. From Table [JJ one 
may observe that it is more likely that our model accommodates the observed dynamics of the universe when 
the number of internal dimensions is higher. As we discussed on Fig. [2] we do not expect a successful model 
for w = 1 and n — 1 and hence these values can be excluded, since g a ,int = 2.46 > 2. Indeed, one may check 
from Table [JJ that in this case the transition to the accelerated expansion is too recent z acc = 0.07 and the 
value of the jerk parameter is too high j a Q = 12.47. In the case for n = 2 too, the transition to the accelerated 
expansion is too recent, although the values q a .mt = 1-74 < 2 and j 0> o = 7.97 are acceptable. We would like 
to discuss the cases n — 6 (l + 3 + n = 10) and n = 22 (l + 3 + n = 26) together. With the choice w = 1 
of this section, the action we considered at the beginning reduces to the gravi-dilaton action of the anomaly 
free superstring theory for n = 6 and of the anomaly free bosonic string theory for n = 22. The deceleration 
parameter of the external space that reaches q a $ = —0.81 at the present age of the universe starts with a 
value q ,int = 1-12 for n = 6 and q a ,mt — 0.85 for n = 22. The jerk parameters, on the other hand, evolves 
from j 0) int = 3.64 at t = to j a fi = 4.97 at itoday = 13.7 Gyr for n = 6 and from j a ^ n t — 2.28 at t = to 
jafl = 3.88 at ttoday = 13.7 Gyr for n = 22. These values can accommodate the observed dynamics of the 
universe [371 GUI EM EH El EH S3- On the other hand, the transition to the accelerated expansion occurs 
more recently z acc = 0.22 in the case for n — 6 to be compared with the value z acc = 0.31 in the case for 
n = 22, with the latter case accommodating the observations better [61 171 [81 [51 ITD1 ITT] . 
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Figure 3: Scale factors of the external and internal spaces 
versus cosmic time t for different number of internal di- 
mensions. We plotted for n = 1,2,6,22 and for n — s> oo 
limit, in the order from the solid blue curve to the dotted 
red curve. 
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Figure 5: Deceleration parameters of the external space 
versus cosmic time t for different number of internal di- 
mensions. We plotted for n = 1,2,6,22 and for n — s> oo 
limit, in the order from the solid blue curve to the dotted 
red curve. 
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Figure 4: Hubble parameters of the external and internal 
spaces versus cosmic time t for different number of internal 
dimensions. We plotted for n = 1,2, 6, 22 and for n — 5- oo 
limit, in the order from the solid blue curve to the dotted 
red curve. 
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Figure 6: Jerk parameters of the external space versus 
cosmic time t for different number of internal dimensions. 
We plotted for n = 1,2, 6, 22 and for n — > oo limit, in the 
order from the solid blue curve to the dotted red curve. 
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Figure 7: Deceleration parameters of the external space versus cosmic redshift z for different number of internal 
dimensions. We plotted for n = 1,2, 6, 22 and for n — > oo limit, in the order from the solid blue curve to the dotted 
red curve. 
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From Table[T]we also observe that some primordial nucleosynthesis occur in all the models, but T ~ 80 keV 
is reached at different times for different number of internal dimensions. It is interesting that in the limit 
n — > oo, the model approaches right to the t < 887 s limit with a value 884 s. In all the cases, the model 
exhibits almost power-law dynamics at energy scales T ~ 80keV, hence we can safely use the condition 
(|46|). which was obtained by considering a primordial nucleosynthesis that yields acceptable light elements 
abundances with q a ~ 0.82. Substituting w = 1 in (|46p we find that when n — 30 the model can represent 
the universe with a great success starting from primordial synthesis times to the present time of the universe. 



5 The model for w = \ (A = 2) 

In this section we give the exact solution for the case w = \ (A = 2), which is valid provided that £ > and 
Uq > 0. For convenience we redefine the time parameter as t — > t — so as to set a = at t — 0. Now, 
using ([55)) in (f2"5|) we obtain the dilaton field as follows: 

1, ( u, 2 



ip = --\i\[--^-t z + ^/2£t). (62) 



The self-interaction potential of the dilaton field 



11=—= — U ° _ . (63) 
-¥/ t 2 + V2£ t 



The scalar curvature of the higher dimensional spacetime is as follows: 



^=07 T2- (64) 

/ t 2 - V2£ tj 

We obtain the scale factors of the external and internal dimensions as follows: 

^ jF? arctanhf -%£L t-l) , - ^ -^SL arctanhf -%L t- 1) , nr , 

a = a 2 e 3 V J and s — s 2 e " / , (65) 

respectively, where and s-i are integration constants that obey a^s-i 2 = Vq. The Hubble parameters of the 
external and internal dimensions are 



H a = =-J^ and H. = -±i[^ (66) 

3V 3 + n -U t 2 + 2V2£t n V 3 + n -U t 2 + 2^28 t V ' 

The deceleration parameters of the external and internal dimensions are 



We also give the jerk parameter for the external dimensions only, since it is important while discussing the 
model from the observational point of view: 



3 = ^ = iy—? + 3VV— ---18)^-3vW^ + - + 13. (68) 
aHl 2 £ n \ V n n I ^/2£ V n n 

The evolution of the universe starts with a zero size external space and infinitely large internal space at 
t = 0. The external space expands until it becomes infinitely large at a finite time t cn( j = 2 "{jf i while the 
internal space contracts to unobservable scales and then hits a singularity at i C nd- The external space enters 
into the accelerated expansion phase at t acc = ^ (^/2 — and into the super-exponential expansion 

phase at t sc = -^p- We note that all the cosmological parameters are determined by the values of and 
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Table 2: Cosmological parameters for the case w — | (A = 2) for different number of internal dimensions n assuming 
9a,o = -0.81 at ttoday = 13.7 Gyr 





n = l 


n = 2 


n = 6 


n = 22 


n — > oo 


^acc 


0.07 


0.12 


0.22 


0.31 


0.34 


13.7 Gyr - t acc 


2.36 Gyr 


3.01 Gyr 


3.95 Gyr 


4.58 Gyr 


4.91 Gyr 


tend 


28.51 Gyr 


28.81 Gyr 


29.25 Gyr 


29.55 Gyr 


29.70 Gyr 


Qa,'mt-> Ja,int 


3.90, 34.30 


2.87, 19.38 


2.00, 10.00 


1.61, 6.80 


1.45, 6.65 


Qa,0: ja,0 


-0.81, 12.48 


-0.81, 7.98 


-0.81, 4.98 


-0.81, 3.89 


-0.81, 3.48 


Qa,cndi Ja,end 


-5.90, 63.70 


-4.87, 42.62 


-4.00, 28.00 


-3.61, 22.47 


-3.45, 20.35 


so 


< 10~ 43 m 


< 10~ 30 m 


< 10~ 21 m 


< nr 18 m 


< nr 16 m 




10" 24 s < 887 s 


KT 16 s < 887 s 


10" 8 s < 887 s 


10" 5 s < 887 s 


10" 3 s < 887 s 



n only; in particular determines the lifetime of the universe, n determines the initial and final values 
of the deceleration parameters and they together determine the time of the onset of the acceleration of the 
external space. Hence, in principle, in this model also we are able to obtain the number of internal dimensions 
n and ^= upon obtaining the values of the deceleration parameter at the two different cosmic redshifts or 
two different times from observations. As in the previous section, choosing q a _o — —0.81 and t = 13.7 Gyr 
considering the latest observations, we are able to depict the dynamics of the universe for different number 
of internal dimensions n. Accordingly, we plot the scale factors of the external and internal dimensions in 
Fig. [5] and the Hubble parameters in Fig. [5] for n = 1, 2, 6, 22 and also for n — >• oo limit. We plot the 
dimensionless deceleration parameter of the external dimensions in Fig. [TO] and the jerk parameter in Fig. [TT1 
The evolution of the deceleration parameter of the external space with respect to the cosmic redshift z, 



■l + Ve 



tanh 



VG /3 + n 



ln(z + 1) + arctanh 



V6 V 3 + 



-(l + q a (z = 0)) 



(69) 



is depicted in Fig. [TOJ We summarize the plots by giving the cosmic redshift (z acc ) and the time passed since 
the time of the onset of acceleration (13.7 Gyr — i acc ), the time of the Big Rip (f cn d) an d the initial (q a ,int, 
Ja,int), present (q a ,0t 3a,o) an d final (q a , C nd, ja,end) values of the deceleration and jerk parameters in Table 
[5] In Table ® assuming that the presence of matter source does not alter the dynamics considerably, we 
also calculated the time £t~so koV to see whether T ~ 80 keV < 887 s is provided to ensure some primordial 
nucleosynthesis. The maximum size of the internal dimensions today sq is also calculated considering the 
condition sns ^ ^proton in Table (5) 

We observe that while the number of the internal dimensions alters the expansion rate of the internal 
dimensions slightly, it alters the contraction rate of the internal dimensions drastically. We note that the 
deceleration parameter of the external space yields a form that is linear in time with a negative slope. This form 
of the deceleration parameter was suggested phenomenologically by Akarsu and Dereli |49, 50] to generalize 
the cosmological models with constant deceleration parameter so as to obtain kinematics that are consistent 
with the observations. In that study, both the initial value and the slope of the deceleration parameter were 
arbitrary. On the other hand, here the initial value of the deceleration parameter is determined solely by the 
number internal dimensions n, namely q a ^ n t = 2y6 — 1, — 1 for n = 0, oo. From Tabled one may 
observe that it is more likely that our model fits the observed dynamics of the universe when the number of 
the internal dimensions is higher. Accordingly, the cases n = 1 and n — 2 can be excluded, since q a ,int > 2. 
One may check from Table [2] that the transition to the accelerated expansion is indeed too recent and the jerk 
parameters evolves with very high values up to the present time of the universe. We would like to note an 
interesting property for the case of n — 6 that corresponds number of spacetime dimensions in anomaly-free 
superstring theories. In this case the initial value of the deceleration parameter of the external dimensions 
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Figure 8: Scale factors of the external and internal spaces 
versus cosmic time t for different number of internal di- 
mensions. We plotted for n = 1,2,6,22 and for n — s> oo 
limit, in the order from the solid blue curve to the dotted 
red curve. 
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Figure 10: Deceleration parameters of the external space 
versus cosmic time t for different number of internal di- 
mensions. We plotted for n = 1,2,6,22 and for n — s> oo 
limit, in the order from the solid blue curve to the dotted 
red curve. 
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Figure 9: Hubble parameters of the external and internal 
spaces versus cosmic time t for different number of internal 
dimensions. We plotted for n = 1,2, 6, 22 and for n — 5- oo 
limit, in the order from the solid blue curve to the dotted 
red curve. 
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Figure 11: Jerk parameters of the external space versus 
cosmic time t for different number of internal dimensions. 
We plotted for n = 1,2, 6, 22 and for n — > oo limit, in the 
order from the solid blue curve to the dotted red curve. 




Figure 12: Deceleration parameters of the external space versus cosmic redshift z for different number of internal 
dimensions. We plotted for n — 1,2, 6, 22 and for n — > oo limit, in the order from the solid blue curve to the dotted 
red curve. 
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gets exact integer values g a! i n t = 2 and j a .int = 10 which correspond to values of the deceleration parameter 
and jerk parameter in the presence of a stiff fluid in general relativity that was proposed by Zeldovich |47) and 
then Barrow [?S] as the most probable EoS state parameter of the very early universe. However, in this case 
the transition redshift is quite recent z acc = 0.22 but it can still be accommodated within observational studies 
(e.g. [HI El [5J [5] ■ In the case for n = 22, both the transition redshift to the accelerated expansion z acc = 0.31 
and the value of the jerk parameter j a ^ — 3.86 are in line with many observational studies [5J [7J [5J [H] • 

In this solution too, the models can accommodate the observed dynamics [61 171 151 151 ITUl ITT1 13"T1 [551 1551 HUl 
[U 021 S3] for the numbers of internal space larger than n = 6. From Table [2] we observe that some primordial 
nucleosynthesis occur in all these models, but T ~ 80 keV is reached at different times for different number of 
internal dimensions. However, we note that the time when the temperature T ~ 80keV is reached too early 
in all the models. Indeed we note that none of the models accommodate the required value of the deceleration 
parameter q a ~ 0.82 at the early times for a primordial nucleosynthesis that yields acceptable light element 
abundances. 



6 The universe according to an observer in (1 + 3)-dimensions 

We argued above that the presence of internal dimensions should not affect the local physical processes, 
for instance the primordial nucleosynthesis, and hence cannot be observed directly and locally, for n > 2 if 
they remain at Planck length scales today and for n > 7 if they remain at just below the scales probed by 
the LHC experiments today. Hence we can consider an observer who lives in the higher dimensional steady 
state universe we constructed in this paper but is unaware of the presence of the internal dimensions. Let 
us now assume that this observer considers the Einstein's general theory of relativity as the valid theory of 
gravity at cosmological scales. Because the observer cannot perceive the internal dimensions she/he would 
represent the geometry of the universe she/he perceives with the (1 + 3)-dimensional part of the metric ([2]) we 
considered, i.e., with the spatially flat Robertson- Walker spacetime. Using the conventional general relativity 
by considering a comoving perfect fluid that would control the dynamics of the observed universe, the observer 
would obtain the following field equations: 

•2 • 2 

3% = kp and % + 2- = -kp, (70) 
a z a z a 

where k = 8nG with G being the gravitational constant scaled in 4-dimensions and p and p are the inferred 
effective energy density and pressure by the observer. In cosmology we can usually characterize a fluid with 
a simple equation of state parameter, which is the ratio between the pressure and energy density of the fluid 
and is not necessarily constant. Hence, the observer is further able to obtain the equation of state parameter 
of the inferred effective fluid using the energy density and pressure: 

w^^ = -\ + \ qa . (71) 
p 3 3 

We apply this argument to the exact solutions we gave for w = 1 and w = 1. 

In the model with £ > and w = 1 (A = 0), using a we obtained for the case w = 1 ([57)) in ([70)) and 

([7l]i. the energy density and EoS parameter of the induced four dimensional effective fluid become 



P=jj— — =- and ™=-l + 2W^ cos( v / %t), (72) 
k (3 + n) snr(V^o t) V 3n 

respectively. The EoS above characterizes a fluid that would give the same kinematics as in general relativity 
for the external space as in our higher dimensional steady state model in dilaton gravity in the source-free 
universe. We plot the energy density of the inferred fluid in Fig. [13] and the EoS parameter in Fig. [14] for 
different numbers of the internal dimensions n by setting q a ^ = —0.81 at itoday = 13.7 Gyr for the present 
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Figure 13: The inferred energy density p by an ob- 
server living in (1 + 3)-dimensions versus cosmic time t 
for different number of internal dimensions. We plotted 
for n — 1, 2, 6, 22 and for n — > oo limit, in the order from 
the solid blue curve to the dotted red curve. 



Figure 14: The EoS parameter of the inferred fluid by 
an observer living in (l + 3)-dimensions versus cosmic time 
t for different number of internal dimensions. We plotted 
for n = 1, 2, 6, 22 and for n — s> oo limit, in the order from 
the solid blue curve to the dotted red curve. 



age of the universe. Using q a ^ — —0.81 in (|7ip we find that the present value of the EoS parameter of the 
inferred effective fluid is wq — —0.87, which is consistent with the observational studies. However, we note 
that the fluid in fact exhibits a quintom behavior that can evolve below the phantom divide line w = — 1; 
namely, w — 1 + atl a — and w -> — 1 — 2^J a s a — >• oo (See Ref. [5TJ [52] for crossing the 

phantom divide line and quintom cosmology). For instance, the EoS parameter evolves from w a =o = 0.22 
to w a =oo = —2.22 in the case n = 30 that is able to give complete picture of the universe with correct light 
element abundances and accommodating the observed dynamics of the universe. 

In the model with £ > and w = \ (A = 2), using a we obtained for the case w — \ (|65|) in (|70p and 

(|7ip , the energy density and EoS parameter of the induced four dimensional effective fluid become 



~ 1 4n£ , . , „ 3 + nU „ /- /3 + n 

5=^ = and w = -l-2\ £=t + 2V2\ , 73) 

fc(3 + n)(-[/ t 2 + 2 v / 2£t) 2 V 3n y/£ V 3n 

respectively. The EoS here characterizes the fluid that would give the same kinematics as in general relativity 
for the external space in our higher dimensional steady state model in dilaton gravity in a source- free universe. 
We plot the energy density of the inferred fluid in Fig. [TS] and the EoS parameter in Fig. [TH] for different 
numbers of the internal dimensions n by setting q a o = —0.81 at ttoday = 13.7 Gyr for the present age of the 
universe. Using q a ^ = —0.81 in (|TTj) we find that the present value of the EoS parameter of the inferred 
effective fluid is wq = —0.87, which is consistent with the observational studies. In this solution we find the 
case n — 6 interesting due to two reasons: (i) the initial value of the EoS parameter is 1 that corresponds to a 
stiff fluid that was already proposed as the most probable EoS parameter of the very early universe [171 HE] > 
(ii) 1 + 3 + n = 10 corresponds to gravi-dilaton action of the anomaly-free effective superstring theories. 

7 Final Remarks 

We considered (1 + 3 + n)-dimensional cosmological solutions in which the size of the (3 + n)-dimensional 
total volume is preserved in dilaton gravity with arbitrary dilaton coupling constant (w > 0) and exponential 
dilaton self-interaction potential in the string frame. We obtain various classes of solutions: 

• Static universe: Both the external and the internal spaces are non-dynamical. 

• Power-law dynamics: Both the external and the internal spaces exhibit power-law behavior. The 
accelerated expansion of the external space is possible depending on the values of w and n. 
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Figure 15: The inferred energy density p by an ob- 
server living in (1 + 3)-dimensions versus cosmic time t 
for different number of internal dimensions. We plotted 
for n — 1, 2, 6, 22 and for n — > oo limit, in the order from 
the solid blue curve to the dotted red curve. 



Figure 16: The EoS parameter of the inferred fluid by 
an observer living in (l + 3)-dimensions versus cosmic time 
t for different number of internal dimensions. We plotted 
for n = 1, 2, 6, 22 and for n — s> oo limit, in the order from 
the solid blue curve to the dotted red curve. 



• Evolving deceleration parameter solutions: Both the external and the internal spaces exhibit time 
varying deceleration parameters. Depending on the values of w and n, there are cases in which the 
external space exhibits late time acceleration. Although it was not possible to obtain analytical solutions 
for arbitrary values of the dilaton coupling constant w, we were able to give a discussion on the general 
behavior of such models since the value of the deceleration parameter when the expansion of the external 
space starts is determined by w and n only. We give explicit and exact solutions for w = 1 and w = h. 
Assuming the present value of the deceleration parameter qo — —0.81, the redshift of transition to 
accelerated expansion is 0.07 < z acc < 0.34 for w = \ and 0.07 < z acc < 0.37 for w = 1 from one 
internal dimensions to infinitely many internal dimensions. These are encouraging results. The best 
w and n combination that accommodates the observational data would be obtained easily if we were 
able to obtain solutions analytically for arbitrary values of w. Nevertheless, using numerical methods 
we still seek for the best w and n combination. Once this is achieved one would be able to give an 
estimate about the number of the internal dimensions and the time for the Big Rip. The Big Rip end of 
the universe is also in accordance with the many studies that argue that the cosmological observations 
favor phantom cosmologies [29j |32l [31] rather then quintessence and ACDM cosmologies. 
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